LOAN COPY: RETURN TO 
AFWL (WLIL-2) 
KIRTLAND AFB, N MEX 


A RENDEZVOUS AND STATION-KEEPING 
SIMULATOR WITH FLY-AROUND CAPABILITY 
AND VIRTUAL-IMAGE VIEWING 


by William M. Kahlbaum, Jr., Richard E. Bardusch, 
and John J. Stumpf 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


APRIL 1969 



TECH LIBRARY KAFB, NM 



□ 1318 


•n 


A RENDEZVOUS AND STATION-KEEPING SIMULATOR 

WITH FLY- AROUND CAPABILITY AND 

VIRTUAL-IMAGE VIEWING 

By William M. Kahlbaum, Jr., Richard E. Bardusch, 
and John J. Stumpf 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — CFSTI price $3.00 



A RENDEZVOUS AND STATION-KEEPING SIMULATOR 
WITH FLY- AROUND CAPABILITY AND 
VIRTUAL- IMAGE VIEWING 

By William M. Kahlbaum, Jr., Richard E. Bardusch, 
and John J. Stumpf 
Langley Research Center 

SUMMARY 

A screen-type-refractive-virtual-image rendezvous simulator has been designed 
and built at the Langley Research Center. In this simulator, the test subject views com- 
bined images of a celestial background, a featureless earth horizon, and a rendezvous tar- 
get vehicle through a circular window with a 30° field of view. The virtual images of the 
horizon and stars appear to be at infinity, while the target-vehicle image appears to be at 
some closer distance. The three elements of the display are rear projected on screens 
by means of separate gimbaled projectors. The horizon and star images are generated 
optically, while the target-vehicle image is generated by a closed-circuit television sys- 
tem. The observed size of the target-vehicle image varies as a function of problem range. 
Range changes within 160 meters of the observer are simulated by moving the television 
camera along a range bed. Range changes beyond 160 meters with satisfactory television 
line resolution are simulated by continuously shrinking the raster on the television projec- 
tion tube. Since the problem requires three degrees of continuous rotational freedom for 
the observer, a central computer program generates attitude information (by means of 
quaternion rate equations) and presents this information to the simulator in the form of 
direction cosine matrices. Special-purpose resolver computers are employed to convert 
the direction cosines into specific hardware gimbal angles. 

This report describes the simulator equipment and gives particular attention to the 
raster shrinkage and gimbal-angle resolving techniques. 

INTRODUCTION 

A simulator has been developed at the Langley Research Center to provide a real- 
istic visual display for use in rendezvous and station-keeping studies. The visual display 
furnishes the observer with a dynamic TT out-the-window TT scene consisting of a target- 
vehicle image, a featureless earth horizon, and a star-field background. Dynamic motion 
of the scene is achieved by means of gimbaled projection devices. A photograph of the 



combined image of the horizon and target as viewed through the window is presented as 
figure 1. Further visual realism is achieved by the use of a screen-type-virtual-image 
viewing system. 

The original intent in the development of this simulator was to investigate the 
optical characteristics of screen-type-virtual-image systems. Shortly after its con- 
struction, however, the system was adapted for use in rendezvous and station-keeping 
studies for the Gemini program. Although the overall system concept is not unique, two 
items of the equipment developed represent techniques which are considered to be of 
general interest and applicability in the field of simulation. The first is a raster shrink- 
age technique used with an 80-kilovolt, theater-type television projection system. The 
second is a technique for employing resolvers to convert direction cosines into gimbal 
position angles. This second technique is particularly attractive in that the resolver 
equipment can be programed easily for any order of gimbal rotation which may be 
required by the projection hardware and can provide continuous rotational freedom about 
all axes. 

SYMBOLS 

A amplifier 

a radius of horizon mask, inches (meters) 

ay elements of 3-by-3 direction cosine matrix 

b distance from center of rotation of gimbal ed mirror to projection screen, 

inches (meters) 

C capacitor 

Cfc capacitance of any element k, microfarads 

D lateral displacement of target image on projection screen relative to system 

optic axis, inches (meters) 

D’ lateral displacement of target virtual image relative to system optic axis, 

inches (meters) 

d distance from television projection lens to center of rotation of two-gimbal 

mirror, inches (meters) 
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voltage appearing on capacitor Cl, volts 


voltage output of amplifier A2 plus a bias voltage proportional to target 
range, volts 

voltage proportional to target range from central computer, volts 

voltage output of amplifier Al, volts 

voltage output of amplifier A2, volts 

voltage output of amplifier A3, volts 

horizon-projector mask height, inches (meters) 

focal length of television camera lens, inches (meters) 

focal length of television projection lens, inches (meters) 

focal length of viewing lens, inches (meters) 

distance between television camera lens and image formed on face of 
vidicon, inches (meters) 

lateral dimension of horizon image on screen relative to system optic axis, 
inches (meters) 

current, amperes 

lateral dimension of real-world target relative to line of sight, feet (meters) 

lateral dimension of target virtual image relative to system optic axis, 
inches (meters) 

lateral dimension of target model relative to camera optic axis, inches 
(meters) 

lateral dimension of image of target model on face of television camera 
pickup tube, inches (meters) 



lateral dimension of target image on face of television projection tube, 
inches (meters) 

lateral dimension of target image on projection screen, inches (meters) 

nominal distance from viewing lens to observer's eye position, inches 
(meters) 

distance from gimbal center of star projector to projection screen, inches 
(meters) 

distance from gimbal center of horizon projector to projection screen, 
inches (meters) 

lens magnification, dimensionless 
electronic magnification factor, dimensionless 

distance from television projection lens to target projection screen, 
inches (meters) 

transistor 

resistor 

resistance of any element k, ohms 
target-to- observer range, feet (meters) 

object distance between television camera lens and center of rotation of 
target- model gimbal, inches (meters) 

distance from projection screen to viewing lens, inches (meters) 

distance from viewing lens to observed virtual image, inches (meters) 

time, seconds 

displacement of a projected star image on screen relative to system optic 
axis, inches (meters) 


a 

a' 

|8 

e 

Oh 

oi 

6 r p 

<h 

0T 

4>o 

^1 

l//rp 


line-of- sight yaw angle, radians 
two-gimbal- mirror yaw drive angle, radians 
line-of-sight pitch angle, radians 
two-gimbal- mirror pitch drive angle, radians 

observed angular displacement of a star relative to system optic axis, 
radians 

horizon-projector pitch gimbal angle, radians 
star-projector pitch gimbal angle, radians 
target- model pitch gimbal angle, radians 

observed horizon depression angle relative to system optic axis, radians 

angular displacement of a star on star projector sphere relative to system 
optic axis, radians 

horizon-projector roll gimbal angle, radians 
star -projector roll gimbal angle, radians 
target-model-image roll gimbal angle, radians 

angle subtended by horizon projector mask relative to system optic axis as 
measured at gimbal center of horizon projector, radians 

horizon- projector yaw gimbal angle, radians 

star-projector yaw gimbal angle, radians 

target- model yaw gimbal angle, radians 

orthogonal unit vectors defining observer body-axis system Xj,Yj,Zj, 
dimensionless 
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orthogonal unit vectors defining line-of-sight axis system X2,Y2,Z 2 , 
dimensionless 

ig, [3,6:3 orthogonal unit vectors defining target body-axis system X3,Y3,Z3, 
dimensionless 

14,34, £4 orthogonal unit vectors defining observer local vertical-axis system 

X4,Y4,Z4, dimensionless 

{5, [5, 65 orthogonal unit vectors defining star-field (inertial) axis system X^YgjZg, 

dimensionless 

H amplifier open- loop gain, dimensionless 

DESCRIPTION OF OVERALL SYSTEM 

The physical system consists of two vehicles (the observer and the target) each of 
which is represented in the computer mathematical model by equations of motion in six 
degrees of freedom. Figure 2 illustrates the angular relationships between the several 
axis systems which describe the physical system. The simulation employs a fixed-base 
cockpit so that all scene motions must be derived relative to the observer body-axis 
system. 

A block diagram of the entire simulator-computer system is presented as figure 3. 
The cockpit is equipped with an instrument panel and controllers with the pilot's con- 
troller deflections being fed to the computer as inputs to the equations of motion. The 
computer outputs generate the signals to drive the cockpit displays and projection sys- 
tems, which form the images of the target, horizon, and star field. As shown in figure 3, 
the line-of-sight angles a and /3 are directly available from the computer program 
as inputs to the target positioning mirror. All other gimbal angles exist implicitly as 
elements of the direction cosine matrices representing the several axis transformations. 
Special-purpose resolver computers (Euler angle computers) are used in the conversion 
of the matrix elements to the required gimbal angles with the angular outputs generated 
in the form of shaft rotations. A synchro device attached to each shaft provides the posi- 
tion input information to the closed-loop servo drive on each projector gimbal. This 
method of generating the gimbal angles provides the fly- around capability which was 
required during the station- keeping portion of the simulation. The target image is 
generated by means of a closed- circuit television (T V ) system. The range excursion 
of 10 meters to 3600 meters necessitated the use of a combination of physical camera 
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motion (range bed) and continuous shrinkage of the projection-tube raster to provide a 
properly scaled image with satisfactory line resolution at all range values. 

VISUAL DISPLAY 


Formation of Virtual Images 

Figure 4 illustrates the projection and optical mixing of the star-field, horizon, 
and target-vehicle images. The images on the rear projection screens SI and S2 are 
combined by means of beam splitter B2 and appear as a composite image to the observer. 
Screen SI, which displays the star-horizon scene, is located at the focal point of the 
viewing lens and, thus, a virtual image is formed at infinity. Screen S2, which displays 
the target image, is located between the lens and its focal point and, thus, a virtual image 
is formed at a predetermined finite distance. In either case, the apparent distance of 
the virtual image may be determined by using the following formula for a simple lens: 


where 


S' = 


Sf. 


v 


S - f„ 


( 1 ) 


S distance from projection screen to viewing lens 

S' distance from viewing lens to observed virtual image 

f v focal length of viewing lens 

Since the apparent distance of the target image is determined by the position of screen S2 
relative to the focal point of the viewing lens, it would be preferable to vary this relation- 
ship as a function of target range in order to maintain the correct observed parallax. 
Although this capability does not exist in the equipment, screen S2 is equipped with a 
mechanism which allows the adjustment of the apparent image distance. 

Star- Field- Image Generation 

General description. - The star field from a point- light- source planetarium pro- 
jector is imaged onto screen SI after reflection from beam splitter Bl. The planetarium 
projector contains a total of 1500 accurately positioned stars. Sixty of the brighter stars 
are projected optically by means of small lenses mounted in the surface of the sphere. 
The remaining stars are produced by small holes drilled in the sphere which allow the 
light from the point source to fall on the screen; star magnitudes are determined by the 
size of the holes. (See fig. 4.) 
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Motion .- The order of rotation of the star-field projector gimbals is <pj, 0j, \f/j. 
As shown in figure 2(c), the orientation of the rotating star-field (inertial) axis system 
relative to the fixed-body axis system is specified by three Euler angle rotations: <pj 

about the i^-axis, followed by 0j about the intermediate j^' -axis, followed by about 
the k 5 -axis. The resulting transformation is (see ref. 1) 


V 


C0jC^j 

C 4> jS + s <p jS 0jC 

S<fijS - C0jS0jC^j 

Y l~ 


Y 5 

= 

- c 6 jS 

C0JC1//J - S0jS0jSl//j 

+ C(f)jS ^jS^j 

Y 1 

(2) 

_ Z 5_ 


s 

- S0JC0J 

C 0jC 

_ Z L 



where c and s represent cosine and sine, respectively. 

Projector geom etry.- The angular relationships between the stars on the planetar- 
ium sphere must be preserved in the observed star -field image. The requirements 
which must be satisfied to guarantee this relationship are illustrated in figure 5. Con- 
sider two stars A and B which are separated by an angle <p on the sphere (star A is 
assumed to be coincident with the optic axis for simplicity). If it is assumed that all 
light rays from a given star emerge from the viewing lens parallel to each other, the 
following expressions can be written: 


tan 4 > = 

L 

(3) 

tan 8 = t - 
tv 

(4) 


where 8 is the observed angle between the two stars. The required condition for 
8 = <fr is that L = f v which specifies the distance of the star projector from the screen. 

Horizon-Image Generation 

General description .- Figure 6 illustrates the geometry involved in the projection 
of the horizon image. A point light source is mounted at the center of the cylindrical 
mask so that a shadow is formed on screen SI. The constant height of the horizon- 
projector mask e determines the depression angle and curvature of the horizon image 
on the screen. Relative motion of the horizon in roll and pitch is introduced by rotating 
the mask in a two-gimbal system. 

Motion .- The order of rotation of the horizon projector gimbals is 0^, 0^. As 

shown in figure 2(d), the orientation of the rotating local vertical-axis system relative to 
the fixed-body axis system is specified by two Euler angle rotations: <p ^ about the 

ij-axis followed by 6^ about the J^j-axis. The third rotation iph about the l^-axis is 
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not included since the horizon image contains no terrain features. The resulting trans- 
formation is 


xf 


ce h 

s0 h se h 

-c0hS0 h 

xf 

Y 4 

= 

0 

C0 h 

s0h 

Yl 

_ Z 4_ 


se h 

-s0 h c0 h 

°0h c0 h_ 

_Zl_ 


(5) 


Projector geometry .- The observed depression angle of the horizon display is 
established by the height of the horizon-projector mask e with the required depression 
angle being determined by the desired altitude. This relationship is illustrated in fig- 
ure 6 for 0^ = 0. If it is assumed that all light rays from a given point on the horizon 
emerge from the viewing lens parallel to each other, the following expressions can be 
written: 

H 0 = f v tan 6 0 = M tan <p Q (6) 

or 

f 

tan <P 0 = jj [ tan 0 o ( 7 ) 

Also, from figure 6 

tan 0Q = | (8) 

Combining equations (7) and (8) gives 

| M = f v tan e Q 

This expression determines the mask height and location for any desired depression 
angle. 


Target-Vehicle-Image Generation and Display 

General description .- The target-vehicle image projected on screen S2 is generated 
by means of a closed-circuit television system as illustrated in figure 7. The television 
pickup consists of the camera lens L3 and a vidicon pickup tube. The vidicon and its 
control unit operate at 675 lines per frame, 30 frames per second with a 2:1 interlace. 

The vidicon and its deflection yoke are mounted inside a servo-driven slip- ring assembly, 
the motion of which introduces the roll of the target image about the line of sight. The 
remaining two motions of the target vehicle relative to the line of sight are produced by 
the two gimbal rotations of the target model. The camera and slip- ring assembly are 
mounted on the movable carriage of a range bed. The size of the image formed on the 
face of the vidicon is proportional to r'. Proper focus is maintained by varying the 
effective focal length of the camera lens as a function of r\ Any value of target range 
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beyond the maximum limit of the range bed is simulated by varying the electronic magni- 
fication factor m e in the projection electronics. The magnification factor is defined as 

m _ Size of projection-tube raster 
e Size of vidicon-tube raster 

Since the effects of varying m e are so closely related to the subject of raster shrink- 
age, detailed discussion of this subject is deferred to the section entitled "Target- 
Vehicle-Range Scaling and Raster Shrinkage." 

The target vehicle image is formed on the face of a high- intensity projection kine- 
scope which is capable of a maximum highlight brightness of 8000 foot-lamberts 
(27 410 candelas/ meter^) when the background raster is extinguished. The target image 
is then projected onto screen S2 by means of lens L2 and the two-gimbal mirror. The 
measured highlight brightness on the screen is approximately 5 foot-lamberts 
(17 candelas/ meter 2) . The two-gimbal mirror positions the target image in azimuth and 
elevation, which correspond to the line-of-sight (LOS) angles a and /3. 

Motion .- The relative rotational motion between the observer body-axis system 
and target-vehicle body-axis system may be broken down into two distinctly separate 
types as illustrated in figures 2(b) and 2(a). First, there is rotation of the LOS axis sys- 
tem relative to the observer body-axis system which is completely specified by the LOS 
angles a and /3. Second, there is rotation of the target body-axis system relative to 
the LOS axis system which is defined by the three Euler angle rotations: cpj, about the 

i2~axis, followed by 0< p about the intermediate j2 ~axis, followed by i//p about the 
kg-axis. (Note that the angles 0-p, i/'t* anc * 0T fully describe target aspect with 

respect to the LOS.) The relative rotational motions that result from basic vehicle 
motions are described by the following examples which are presented as an aid in visual- 
izing these motions: 

•v 

(1) Relative translation between the two vehicles results in variations of a, j3, 

0-p, d'Y, and i^rp. However, for certain relative orientations, <p ' p does not vary (for 
example, target and observer body axes parallel). 

(2) Rotation of the observer vehicle results in variation of the angles a, /S, and 

<£rp. However, target aspect angles (ip^ and 6< p) do not change since the axis system 
i2,j2 does not rotate relative to the target body-axis system i3 , 5^3 , kg . 

(3) Rotation of the target body-axis system results in variations of the angles <p r p, 
£?rp, and ^/ /r p* 

The line-of-sight angles are directly available from the main computer program so 
that no axis transformation is required. The transformation defining the target gimbal 
angles is 
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x 3~ 


C 0rjiC^/rp 

c0t s ^t *^^T^^T^^T 

Stjb'pS^fp - C<prj\S 0rjiCl//rp 

X 2 

Y 3 


- C 0rpSl//rp 

C^^pC^rp - S0rpS0rpSl^/p 

S <^)rp Cl^/rj! + C Qrj*Sy\£/rji 

y 2 

Z3 


s 6 

-S(^>rpC 0rp 

C (j)rpC Qrj* 

Z2 


Projector geometry .- Certain conditions must be satisfied in order to have the 
mirror drive angles a' and j3’ equal to the line-of- sight angles a and /3. From 
figure 8 the following relationships may be written: 


tan a = 


D* 

K + S' 


( 10 ) 


tan « T = ^ 
b 


(ID 


where D f is the lateral displacement of target virtual image relative to system optic 
axis and D is the lateral displacement of target image on projection screen relative to 
system optic axis. From simple lens theory (ref. 2) the lateral magnification of the lens 
may be written as 


m = 


D T 

D 


£T 

S 


( 12 ) 


Since S T » K, equation (10) may be approximated by 


tan a ~ 


D t 

S’ 


D 

S 


(13) 


Comparison of equations (13) and (11) indicates that if a is to equal oi\ then 

b = S (14) 


which fixes the position of the mirror gimbal center relative to the screen. An identical 
derivation using /3 would show that the same requirement holds in order to have /3 = /3\ 


Target-Vehicle-Image Generation Equipment 

Figure 9 illustrates the mechanical arrangement of the range-bed system and, in 
particular, the manner in which the camera assembly is mounted on the range-bed rails. 
A closed-loop servo system is used to position the camera assembly on the rails. Fig- 
ure 10 shows a close-up view of the camera assembly, range-bed rails, and the target 
model mounted in the two-gimbal system. Figure 11 illustrates the relative locations of 
the TV projection tube, projection lens, and two-gimbal mirror. The projection- tube 
housing is lined with lead for protection from the X-rays generated by the 80-kilovolt 
projection tube. Figure 12 shows an overall view of the entire projection assembly with 
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the rack on the left housing the raster shrinkage equipment. The projector mounting 
frame is equipped with rollers for ease of movement and jack screws for vertical 
adjustment. 


TARGET-VEHICLE RANGE SCALING AND RASTER SHRINKAGE 


Derivation of Range Scaling Equations 

In order to provide the full problem range excursion and still maintain correct 
target perspective at near range, the target- vehicle- image generation system employs a 
combination of physical camera motion (range bed) and projector raster shrinkage. 

Correct perspective requires that the range-bed scale factor be identical with that 
of the target model. Therefore, the ^ -scale model fixes the maximum simulated range 
of the range bed at 160 meters. With this scale factor, the near limit of the range bed 
fixes the minimum simulated range at 10 meters. A technique has been developed which 
extends the maximum range capability of the system while preserving the TV image- line 
resolution that exists at the 160- meter range (maximum range-bed travel). This tech- 
nique consists of continuously shrinking the projection-tube raster from its maximum 
size whenever the problem range exceeds 160 meters. 

The following equation development (see fig. 13) serves to establish a relationship 
between model dimension and observed image dimension for the combined range bed and 
raster shrinkage system. The overall magnification of the television system may be 
determined by repeated application of the simple lens magnification equation (ref. 2) to 
the three lenses in figure 13, as follows: 



Combining equations (15) to (18) results in the following equation which relates the dimen- 
sion of the final image to the corresponding dimension of the target model: 


J’ = 



(19) 
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The dimension J' is related to the corresponding dimension of the real-world object J 
by 

S’ + K = r 

Substitution of equation (20) into equation (19) results in the following expression: 

f 


J(S’ + K) _ / S’ - f v \ /n - f p \ 


( m e) 


r’ - f. 


Jc 


( 21 ) 


Equation (21) may be rewritten to illustrate the relationship between the range-bed range 
and real-world range which is used, with constant m e , for ranges less than or equal to 
160 meters. The resulting equation is 


r' = 



( 22 ) 


where J C/ /J is simply the model scale factor. In equation (22) the term f c is nor- 
mally small in comparison to the other term and may be neglected. Equation (21) may 
also be rewritten to show the relationship between the electronic magnification and real- 
world range which is used, with r' constant, for ranges greater than 160 meters. The 
resulting equation is 




J(S’ + K) 

Jc IS 



(23) 


General Circuit Considerations 

Closed-loop computational elements are utilized to generate the desired voltage 
and current waveforms in a nonresonant deflection amplifier in order to provide the elec- 
tronic means of symmetrically varying the size of a television raster to simulate image 
ranging. Figure 14 is a block diagram of a single channel of the reduction raster deflec- 
tion system. A composite (horizontal and vertical) raster generator requires two similar 
channels. Although the deflection systems exhibit conceptual similarities in the vertical 
and horizontal channels, the power levels, timing, and internal gains are different. With 
the exception of the synchronization circuitry, the function blocks contain closed-loop 
chopper-stabilized computational amplifiers. 

The video and synchronization voltages from the camera control are separate inputs 
to the projector electronics. The projector electronics consist of a high-gain video 
amplifier and a specially designed sweep-drive unit. The sweep-drive unit is designed 
so that the size of the projection-tube raster may be varied continuously under control 
of the main program range signal. Since only the size of the raster is being varied, the 
number of raster lines contained in the image is constant. As the raster size is reduced, 
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the tube beam current must be reduced proportionately to maintain constant image bright- 
ness. This action also prevents damage to the tube phosphor, especially at small raster 
sizes. 


Deflection Circuitry 

The deflection amplifier is a classical feedback circuit which converts the sum of 
the input voltage waveforms to a proportional yoke current waveform. The conversion 
accuracy is dependent upon the deflection driver's open-loop gain-bandwidth response. 
The voltage excursion capability is governed by the rapid retrace speed of the sawtooth 
current passing through the inductive yoke element. 

Full raster size corresponds to a half- axis deflection angle of 11° with deflection 
angle being linearly proportional to deflection current. Image brightness considerations 
require a cathode- ray-tube ultor voltage of 80 kilovolts. The horizontal deflection cur- 
rent for a 10- microhenry deflection coil is 10 amperes (full-size raster). 

The deflection amplifier is driven by two signals: a voltage sawtooth and a negative 
dc bias both of which are proportional to the range input voltage. A yoke current wave- 
form, symmetrical about a zero reference, will be generated if the magnitude of the bias 
offset is equal to one-half the magnitude of the sawtooth voltage. Variation of the range 
input voltage will cause a symmetrical variation of the yoke current, and hence of the 
raster, about a zero-current center. In addition to satisfying the feedback requirements 
of the deflection amplifier, the voltage waveform obtained at the current sensing resistor 
provides the intensity control circuitry with a reference voltage that is directly propor- 
tional to the raster size. 

The analog range input and the synchronization pulses are buffered for noise. A 
common mode rejection circuit is used for the analog input and contains an inverter- 
amplifier configuration which subtracts the common mode voltage from the analog com- 
puter input. In addition, this circuit provides low impedance drive ahead of the voltage 
stages of the raster generator. The signal conditioning circuitry provides a 2. 5- volt 
noise rejection threshold for incoming synchronization pulses. 

Figure 15 illustrates the circuit configuration of the high-speed, resettable, boot- 
strap integrator used to generate the voltage sawtooth sweeps. The integrator is formed 
by a feedback arrangement of two inverting amplifiers. Linear voltage ramps are gen- 
erated by a constant- current charging capacitor Cl. Transistor Q1 operates as a switch 
and provides the discharge current path for capacitor Cl during retrace time intervals. 
Synchronization pulses applied to the base of Q1 define sweep and reset times. 

The following derivation assumes A1 and A2 to be ideal computational amplifiers. 

A current summation at node A for Q1 in a nonconducting state yields 


14 


or 


12 = 13 - h 


With 


. 2 = ^ 


Re 


E c ) E c 
" Ri 


(24) 


“R5E1 RsEr 

Ea = rT 

and 


equation (24) takes the form 


E c 


RfEi 

R 2 


h 


R - 5 - E — - E r + Ej + — = - R5 E r + E-. U 6 -*- — 

R 6 R 3 r 4 r 6 r 2 r 6 r 2 r 4 r 6 1 \ R 2 R 6 



(25) 


r 5 _ r 6 + R 1 

R3 ~ r 2 


equation (25) reduces to 


T = r 5 
2 " r 4 r 6 


E r 


The current I2 may also be expressed as 


I 2 = Ci 


dEc 

dt 


CiRi dE x 

R2 "dT 


Substituting equation (28) into equation (27) results in 


dEj R 2 R 5Er 
dt RjR^RgC j 

or 


_ R 2 R 5 C 

0 R^^gC^ J 


E r dt 


(26) 


(27) 


(28) 


(29) 


Equation (29) defines the output Ej to be a linear voltage ramp for the condition 
of a dc input reference voltage. Ramp linearity is dependent upon the practical realiza- 
tion of resistive equality as defined by equation (26), inequality resulting in either a bow 
or sag in ramp voltage. The most stringent requirement which the practical circuit 
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places upon the ideal configuration is that of bandwidth. The problem is minimized, 
however, in that the amplifiers operate with a near unity-gain feedback loop. The reset 
speed achieved by the horizontal sawtooth generator is of the order of 1 microsecond. 
The feedback integrator afforded a method of referencing the chopper transistor Q1 to a 
circuit ground rather than the virtual ground of an operational amplifier summing 
junction. 

The transfer characteristics of the deflection amplifier are derived for the ampli- 
fier configuration of figure 16. The derivation assumes one input voltage. Node B is a 
summing junction and performs the arithmetic sums of multiple inputs. This fact is 
later used to establish the operating point of the sweep amplifier. 


The input signal is amplified and applied to the yoke and sensing resistor combina- 
tion (L and R 9 ). The yoke current waveform obtained at node C is returned to node B 
by the resistance path R 8 . The subtraction process reduces to zero the arithmetic dif- 
ference between the applied current waveform and the yoke current waveform. 

If the input impedance of amplifier A3 is high and if Rg « Rg, then the voltage 
present at node B can be expressed as 


_ ( E in - l4 R 9) R 8 , R 

E B ~ r 7 + Rg + j 4 r £ 


(30) 


The open- loop gain of amplifier A3 is defined to be -p; therefore, 


% 


( E in " ^4 r 9) r 8 


R 7 + Rg 


+ I4R9 


(31) 


If the open- loop gain of the amplifier is large, then from equation (31) 
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** 8 % \ 
r 7 + r 8 / 


~ R 8 E in 
Rry + Rg 


Reducing this equation gives 

I = ~ R 8 E in 

4 r 7 r 9 


(32) 


Equation (32) defines a linear relationship between the output current and the input 
voltage if the gain of the amplifier is large. Equations (29) and (32), in conjunction with 
yoke test data, are used to establish the voltage and current deflection parameters. The 
gain of the bias amplifier (fig. 17) establishes a current offset equal, to one-half of the 
sweep current. 
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Figure 17 illustrates the basic dc circuit for the horizontal and vertical deflec- 
tion amplifiers. The power stage, shown in transistor form, contains a number of cur- 
rent stages which are represented by Q3 and Q4. Resistors RIO and Rll compensate 
for unequal current gains in addition to providing circuit stability. The amplifier is 
biased by diode D1 to permit a controlled amount of current to pass from transistor Q3 
to Q4. This action minimizes the crossover distortion that is present at ground level. 
Resistor R8 defines the feedback path for both the voltage and current waveforms obtained 
at R9. The power amplifier is driven by a high bandwidth voltage amplifier denoted by 
A3. 

The dc voltage gain of the power amplifier approaches that of a complementary 
emitter follower. The transient response, however, is somewhat less than a comple- 
mentary pair owing to the increased capacitance introduced by the large surface area of 
the heat sinks which were directly connected to the collector to achieve the necessary 
heat dissipation. 

In order to maintain a constant picture brightness, with respect to a varying raster 
size, it is necessary to modulate the intensity of the projection tube with an appropriate 
variable bias voltage. Figure 18 depicts the intensity control circuitry. The vertical 
and horizontal sweep currents are sensed and the summation of these terms is func- 
tionally coupled to the control grid of the cathode- ray tube. It is necessary to provide a 
vertical retrace blanking period to overcome the transient capability of the vertical 
amplifiers; that is, less than 50 microseconds of the available 0.5 millisecond are needed 
for current retrace. By serrating this dc bias voltage with the vertical retrace pulse, 
it is possible not only to provide the required vertical blanking but also to allow capaci- 
tive coupling of the resulting waveform into the control grid of the cathode- ray tube. 

Diode dc restoration inhibits the cumulative waveform sag. Amplifier A4 provides the 
necessary voltage gain for the unblank pulse, and the threshold supply determines the 
initial projection- tube bias. 


EULER ANGLE COMPUTERS 

The gimbal angles for the different projectors are provided to the simulator in the 
form of direction cosine matrices, the matrix elements being generated by means of the 
quaternion rate equations (ref. 3). In this section is discussed the mechanization of the 
special-purpose computers which convert the matrix elements into the required shaft 
angles. 

The transformation matrix for each set of gimbals has been defined previously in 
this report. For all gimbals the order of rotation is roll <£, pitch 6, and yaw ip. The 
general form of this matrix is repeated to illustrate the trigonometric dependence of 
each matrix element upon the gimbal angles as follows: 
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X 


(c0jc i/zjj (cepiS \pi + S(pis6iC\l/i) (scp^s - c^s^ci^) 

X 


Y 

= 

(~c6iS\}/i) (ccpicij/i - s^js^si^i) (s$ 1 ci/'i + c0 1 s6>is^/ 1 ) 

Y 

(33) 

_Z_ 

rotating 

s8\ -S0JC01 C0]C 8\ 

Z_ 

fixed 


where subscript 1 refers to the resolver input and subscript 2 refers to the actual 
resolver shaft position in figure 19. For an arbitrary relative orientation between the 
axis systems, the transformation matrix of equation (33) takes on the following form of 
a matrix of nine direction cosines: 


x" 


Y 

= 

_Z_ 

rotating 


a ll 

a 12 

a 21 

a 22 

a 31 

a 32 



(34) 


where the ay's represent the direction cosines between respective axes in the two sys- 
tems. In particular, the direction cosines which define the axis transformations referred 
to in figure 3 are specified by the elements bjj, cjj, and d^. In equations (33) and (34) 
the same convention of transformation from a fixed axis system to a rotating axis system 
has been maintained. 

If the definition of the fixed and rotating axis systems is identical for the simulator 
and computer mathematical model, then the a^'s in equation (34) are identical with the 
respective matrix elements in the mathematical model. If, however, the mathematical 
model definition of fixed and rotating axis systems is opposite to that of the simulator, 
then the a^'s in equation (34) must be set equal to the elements of the inverse (trans- 
pose) mathematical model matrix. 

Two noteworthy characteristics of the matrix elements are as follows: 

(1) The value of each aij is always between +1.0 and -1.0 which simplifies ampli- 
tude scaling of the analog computer. 

(2) There is no specific order of rotation implied in equation (34) which means that 
the same input matrix may be used to represent simultaneously any or all six possible 
orders of rotation for a three-gimbal system. 

Figure 19 shows the mechanization of an Euler angle computer for the order of 
rotation Each level of computation labeled "A" to "E" is made up of a polar 

resolver (ref. 4) utilizing the mechanization convention illustrated in the lower right-hand 
corner of the figure. Equation (33) shows that the elements of row three of the matrix 
contain sufficient information to determine the two angles <pi and 8\. However, it 
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must be noted that the angle <p j is not properly defined at = 90°. This situation 
corresponds to gimbal lock and is described in a subsequent paragraph. 

The elements ag 2 and a 33 are applied as inputs to the A level of computa- 
tion; this results in the two outputs shown in figure 19. The output of line 1 provides the 
null signal [cos #i(0)j to drive the <p2 servo. Output 2 consists of cos 6\ which com- 
bined with a 3 j = sin 9\ forms the inputs to the B level of computation. The outputs of 
level B are defined by the convention used for level A. To compute the yaw angle i//j, 
the terms cos and sin i/'i must be isolated from among the remaining six terms 
of the matrix. Although elements ajj and a^ seem to be the most obvious choice, 
use of these two elements would introduce an undefined point in identical to that 
which occurs in the <p 2 computation. An alternate solution to the problem is to per- 
form the algebraic operations on the second row indicated in the D and E levels of com- 
putation. (The first row of equation (33) would also satisfy this requirement and the 
choice of either one is arbitrary.) The result of these operations provides the two inputs 
(sin and cos to the C level of computation, the operation of which is identical 
to the operation of the A and B levels. Sign inversions of signals are used where indi- 
cated so that the same convention may be utilized in all levels of computation. 

Even though the matrix elements exhibit no discontinuities, the Euler angle com- 
puter does encounter an indeterminate point at 61 = ±90°. This condition corresponds 
to the normal gimbal lock condition of a three-gimbal system and may be seen in fig- 
ure 19 by observing that the null signal for the (p 2 servo is multiplied by cos 9\. Thus, 
as 61 approaches ±90°, the voltage gradient of the null signal decreases. This is the 
normal situation for a polar resolver and the common remedy is to vary the servo gain 
inversely with cos 9. However, as the 6 = 90° point is reached, the automatic gain 
control is no longer effective in maintaining servo position. At this point the <fi 2 shaft 
may drift away from its last known value. The drift of the <£2 shaft is compensated for 
by an equal and opposite motion of the '^2 shaft since 1P2 and d>2 are coupled 
together through the C and D levels of computation. Any drift which is introduced as 
herein described is automatically eliminated as soon as the system emerges from gimbal 
lock because the $2 shaft once again is providing the correct solution. 

The mechanization described provides a solution which is piecewise continuous. 
When the singularity is encountered, either of two solutions may result upon emerging 
from the singularity. Therefore, if the system must pass in the vicinity of 9 = ±90°, to 
provide predictable behavior it becomes desirable to include logic which restricts the 
solution to only one of the two possible solutions. 
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CONCLUDING REMARKS 


A screen-type- refractive-virtual-image simulator has been designed and built at 
the Langley Research Center and has been applied to space rendezvous simulation for the 
Gemini program. Its applications have been extended to research in general visual tech- 
niques which will be pertinent to future space-flight activities. 

A brief description has been presented of the more prominent features of the simu- 
lator equipment, two of which are worthy of mention. First, a method of continuously 
shrinking the raster on an 80-kilovolt kinescope projector provides the simulator with an 
extended range capability and still maintains target-vehicle- image line resolution. 

Second, the use of direction cosines to define the axis transformations and of Euler angle 
computers to derive the gimbal angles from elements of these matrices has provided the 
system with complete fly-around capability. This simulator has received a high degree 
of pilot acceptance which is due largely to the virtual- image display and the fly-around 
capability. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 18, 1968, 

125-19-06-01-23. 
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Figure L- Combined target and horizon images. L-66-3296 
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Figure 3.- Block diagram of simulator system. 
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Figure 5.- Angular relationships between observer and projector when scene is viewed as a virtual image at infinity. 
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Figure 9.- Range bed. L-66-1636.1 
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Figure 10.- Gimbaled target model and camera assembly. 
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Figure 11.- TV projector showing the projection lens and two-gimbal mirror. 
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Figure 12.- TV projector and raster control unit. 
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Figure 14.- Function block diagram of the projector electronics. 
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